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ABSTRACT 

Emission spectra from magnetars in the soft X-ray band likely contain a ther- 
mal component emerging directly from the neutron star surface. However, the 
lack of observed absorption-like features in quiescent spectra makes it diffi- 
cult to directly constrain physical properties of the atmosphere. We argue 
that future X-ray polarization measurements represent a promising technique 
for directly constraining the magnetar magnetic field strength and geome- 
try. We construct models of the observed polarization signal from a finite 
surface hotspot, using the latest NS atmosphere models for magnetic fields 
B = 4 x 10 13 — 5 x 10 14 G. Our calculations are strongly dependent on the 
NS magnetic field strength and geometry and are more weakly dependent on 
the NS equation of state and atmosphere composition. We discuss how the 
complementary dependencies of phase-resolved spectroscopy and polarimetry 
might resolve degeneracies that currently hamper the determination of mag- 
netar physical parameters using thermal models. 

Key words: polarization - magnetic fields - radiative transfer - stars: at- 
mospheres - stars: magnetic fields - stars: neutron - X-rays: stars. 



1 INTRODUCTION 

Thermal models of emission from neutron star (NS) surfaces and observations of NS cooling 
can be used to constrain the star equation of state (EOS), properties of matter at super- 
nuclear densitites, and phy sics at extreme magnetic field strengths (e.g. 



Prakash et al 



2001 



Yakovlev Sz Pethick 



20041 ). Thermal radiation has been detected from several classes of iso- 



lated NSs, including Anomolous X-ray Pulsars ( AXPs) , Soft Gamma-ray R epeaters (SGRs), 



X-ray Dim Isolated Neutron Stars (XDINS; e.g.. 
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Haberl 



20071 ). pulsars (e.g.. 



Becker &: Pavlov 
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2002 



Marshall fc Schulz 



nova remnants (e.g., 



2002 



De Luca et al 



Sanwal et al 



2002 



2005), and central compact objects in super 



Mori et al. 



2005 



Kargaltsev et al.l 120051 ) . In this 



paper, we focus on emission from AXPs and SGRs, which comprise the population of mag- 
netars, NSs that feature quiescent and bursting emission powered b y the decay of strong 
interi or magnetic fields with surface strengths B ~ 10 14 — 10 15 G (e.g.. IWoods fc Thompson 



2006 



ferred ages, and the lack o 



emission (e.g., 



Kaspi 



and the references therein). There is a near consensus in the literature that the mag- 
netar hypothesis is correct. Evidence favoring the magnetar model includes the energetics, 
flux, and timing properties of SGR and AXP flares, long magnetar rotation periods and in- 
sufficient NS rotational energy to power the observed quiescent 
20041 ). This phenomenology is consistent with internal NS magnetic 
fields strengths of B > 10 15 G. However, other than the large magnetic dipole components 
calculated from timing measurements, there are no direct empirical constraints on external 
magnetar-strength magnetic fields. This is in part due to the lack of features in magnetar 
thermal spectra, a characteristic that is difficult to explain, since broad features have been 
observed in the spectra of NSs with lower magnetic fields and temperatures (e.g., XDINS 
and pulsars). 

Moreover, differences between the observed characteristics of magnetars and other classes 
of high field NSs are still poorly understood. The recently discovered population of high- 
field radio pulsars is similar in most respects to canonical pulsars (they are powered by 
rotational kinetic energy and emit persistently in the radio band). However, the high-field 



radio pulsars have 



Camilo et al- 



ong periods and inferred magnetic field strengths of order B > 4 x 10 13 G 



2000 : 



Gonzalez et al. 



2004 



Kaspi &: McLaughlin 



2005 



Livingstone et al 



20061 ). Several authors have speculated that the distinct behaviors of magnetars and high field 
radio pulsars suggest differences in their magne tic field geometries an d possible evolutionary 



Gonzalez et al. 



20041 ). 



connections between the two populations (e.g., 

One of the keys to understanding the nature of magnetars and their relationship to 
other NS populations is in the production of plausible theoretical models for interpreting 
observations. While significant progress has been made in modeling the spectra of thermal 
radiation from NS atmospheres, the use of such models to explain observational data is still in 
its infancy. In practice, typical soft X-ray magnetar spectra are fit equally well by blackbody 
or atmosphere models plus a power-law. The incorporation of vacuum polarization effects 
into NS atmosphere calculations has improved model fits to observed spectra, and has offered 
a possible explanation for the absence of magnetar spectral lines. In principle, predictions of 
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X-ray Polarization in Magnetar Emission 3 

the beaming pattern of radiation along the surface magnetic field can be coupled to phase- 
resolved spectroscopy to discriminate between NS atmospheres and other emission models. 
However, the relatively low photon count rates from many magnetar sources, as well as the 
degeneracies between the NS magnetic configuration, viewing geometry, and compactness 
make it difficult to extract physical parameters from observations. 

As an alternative, several recent works discuss the evolution of photon polarization states 
in NS magnetospheres, showing that, in the magnetar case, significant linear polarization 



fractions are expected . 



(IHevl fc Shaviv 



Wang k, Lai 



2000 



Measurements of significant X-ray polariz ation, at 2.6 keV and 5 .2 keV, 



possi 



2002 



jly containing a 



Heyl et al. 



2003 



unique signature of the strong magnetic field 



Lai fc Ho 



2003a 



van Adelsberg &: Lai 



2006; 



were performed for the Crab nebula using the OS 0-8 satellite (IWeisskopf et al. 



meas urements confirmed an earlier detection by a sounding rocket experiment (INovick et al 



1976). These 



19721 ). However, as of this writing, no subsequent polarization measurements have been made 



for any object at energies E ~ 0.1 — 10 keV (relevant for thermal magnetar emission). Recent 
advances in instrumentation have stimulated interest in future missions to perform polarime- 



Costa et al. 


2001; 


Kallman 


2004; 


Costa et al. 


2006) 



In this paper, we explore the future role of X-ray polarimetry as a complement to spec- 
troscopy in interpreting observational spectra. We will argue that the combination of po- 
larimetry and spectroscopy can constrain several critical NS parameters, including the tem- 
perature, magnetic field strength and geometry, size of emission region, and mass-to-radius 
ratio of highly mag netized NSs with B ~ 10 12 — 10 15 G. We use the lates t magnetar at 
mosph e re models o f Ivan Adelsberg fc Lail (120061) . and expand on the wo r k of 



Heyl fc Shaviv 



(120021 ) . IHevl et all (120031 ). lLai fc Hoi ([2003a]), and Ivan Adelsberg fc Lail (120061 ) . to compute 
the phase-resolved, observed Stokes parameters from magnetars. We assume that these NSs 
have dipole magnetic field strengths of B = 4 x 10 13 G - 5 x 10 14 G and emit from a re- 
gion centered around the star polar cap with modest opening angle. We confirm that the 
polarization signal from a finite region on a magnetar surface retains important information 
about the strength of the magnetic field, as reported in previous works. We show that this 
signal has a strong dependence on the magnetic field and viewing geometry, and a weaker 
dependence on the NS EOS, emission region size, and atmosphere composition. Finally, 
we argue that polarization measurements can break the degeneracy inherent in inferring 
physical parameters from spectroscopic measurements alone. 

© 0000 RAS, MNRAS 000, 000-000 



4 van Adelsberg & Perna 

Our paper is organized as follows: in $2], we discuss our assumptions and the atmosphere 
models used to calculate the emitted polarization fraction from the NS surface. In §3j we 
describe our methods for calculating the observed polarization signal from an extended NS 
polar cap, including relativistic effects. In §HJ we show the results of our calculations for 
several representative cases. Lastly, in §0 we analyze our results, and compare and contrast 
the information obtained from polarization studies to that from phase-resolved spectroscopy. 



2 PHYSICS INPUTS 



In the canonical pulsar model, the magnetospheres of highly magnetized NSs contain tenuous 
plasma whose density is approximated by the Goldreich- Julian formula, ncj = fB/(ce) ~ 
6.9 x 10 12 /i5i4 cm" 3 , where h = f/(l Hz) and Bu = B/jlO 1 * G) a re the NS rotation fre- 
quency and magnetic field strength (e.g.. IShapiro &: Teukolskylll986l ). The magnetar model 
contains a more complicated magnetic field structure compared to the standard dipole case. 
The presence of twisted magnetic fields (e.g., due to crustal motion and magnetic reconfigura- 
tion) leads to significant magnetospheric currents, which can, in prin ciple, maintain electron - 
positron densities that are orders of magnitude larger than ncj (e.g.. Thompson et al.ll2002l ). 
This suggests the possibility that resonant cyclotron scattering in the magnetosphere distorts 
some of the radiation emerging from the NS surface, producing a NT component at energies 
E ~ 1 — 10 keV. Modeling this emission requires sophisticated scattering calculations that 
depend on the magnetic geometry and magnetosphere structure; significant progress has 



been made in several recent works ( Lyutikov fc Gavriil 



2006 



ReaetaL 



2008|). 



Fernandez fc Thompson 



2007 



In this paper, we focus on thermal photons that emerge directly from the NS surface, 
undistorted by scattering processes in the magnetosphere. In addition, we make the simpli- 
fying assumption that the NS has a pure dipole magnetic field. While the actual magnetar 
field structure is likely to be more complicated, it is unlikely to distort the polarization sig- 
nal, which depends on the field structure at distances much greater than the star radius; far 
from the star, the field is dominated by the dipole component. If the emitted polarization 
is instead determined exclusively by processes at the NS surface or scattering, the vacuum 
polarization signature is destroyed. 

In the work below, we assume that X-ray photons are emitted from a hot region with 
T e Q ~ 5 x 10 6 K, centered around the magnetic pole. The bulk emission from the rest of the 
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X-ray Polarization in Magnetar Emission 5 

star is taken to be at a lower temperature and to contribute negligibly to the observed signal. 
We also assume that the size of the X-ray emission region is much smaller than the NS radius, 
R, with an approximately constant magnetic field normal to the NS surface. The opening 
angular radius of the polar cap is defined to be (3 < 2tt. This geometry is consistent with 
what is observed in several magnetar sources. In quiescence and post-outburst, magnetars 
emit a thermal, soft X-ray component, with effective temperatures kBT c g ~ 0.3 — 0.5 keV 
(along with a non-thermal, hard X-ray power-law). Single and two-component blackbody 



fits to observed, thermal X-ray spectra in 
radiation emerges from a small 



er emission radii Rhh 1 km, imply ing that the 



Tiengo et al. 



2005 



Israel et al. 



20071 ). 



Timing studies of AXPs and SGRs yield rotatio n periods of P ~ 10 s , and imply dipole 



magnetic field strengths B ~ 10 — 10 G (e.g., 



2005 



Camilo et al. 



2007 



Israel et al. 



2007j)E Ivan Adelsberg fc Lail (120061 ) showed that, for 



Kulkarni et al 



2003 



Mereghetti et al 



P > 0.01 s, calculations of the phase-resolved photon Stokes parameters are independent of 
rotation period. 

In §2.11 we briefly discuss radiative transfer in magnetar atmospheres. In §2.21 we give a 
short review of the physics of vacuum polarization, and its effect on magnetar polarization. 



2.1 Magnetar Atmosphere Models 



Emission from magnetar surf aces is de 



; ermined by radiative trans 



Ho et al. 



2003 



van Adelsberg et al. 



er through an atmosphere 



2005 



Medin fc Lai 



2006, 



layer of ionized plasma (see 
for cases in which these conditions may be violated). Due to gravitational setting, this 
atmosphere is likely to be composed of light elements, predominantly hydrogen or helium. 

In a NS atmosphere, photons propagate in two polarization modes that have distinct 
interactions with the plasma medium. Ordinary (O) mode photons are linearly polarized, 
mostly in the plane formed by the photon propagation and magnetic field directions. They 
have an absorption opacity approximately equal to that of photons in non-magnetic plasma. 
Extraordinary (X) mode photons are linearly polarized mostly perpendicular to the plane 
formed by the propagation and magnetic field directions. They have an absorption opacity 



reduced by a factor (Es e /E) 2 , whe re Ebp. 
and E is the photon energy (e.g., 



Meszaros 



158£?i4 keV is the electron cyclotron energy, 



19921 ). The reduction factor in the X mode 



For a complete, current reference list of magnetar observations and measured quantities, see 
http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 
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6 van Adelsberg & Perna 

absorption opacity is due to the suppression of particle motion transverse to the strong 
magnetic field, and is as large as 10 6 at photon energies E ~ 1 keV. X mode photons therefore 
decouple from deeper layers of the atmosphere than O mode photons, and, under typical 
conditions, dominate the emergent radiation. Therefore, magnetar atmos pheres emit strongly 



polarized radia t ion, with 



Ho fc Lai 



2003 



Lai fc Ho 



inear p olarization fractions up to 100% (e.g. 



Pavlov et al. 



1994 



2003bl ). In addition, atmospheric emission is beamed along the 



magnetic field, though the beaming pattern broadens w hen the effects of vacuum po larization 



in strong magnetic fields are taken into account (see Ivan Adelsberg &: Lai 
references therein). 



2006 



, and the 



The atmosphere models of 



van Adelsberg Lail (120061 ) are the first to quantitatively 



incorporate the physics of vacuum polarization in strong magnetic fields, and produce ac- 
curate results for fully ionized atmospheres with field strengths B = 10 12 — 10 15 G. We use 
these radiative transfer calculations to construct the linear polarization fraction, defined: 
nem F°(6)-F*(6) 

llE FO(5) + F£(5y [L) 
where F^(5) and F§{$) are the observed intensities for X and O mode photons emerging 
from a hotspot on the star surface, and 5 is the angle between the photon propagation 
direction and surface normal (see §3.21) . 



2.2 Vacuum Polarization 



When the NS magnetic field strength is greater than the quantum critical field Bq ~4x 
10 13 G (obtained by setting the electron cyclotron energy equal to its rest energy), there are 
significant vacuum contributions to the dielectric tensor and magnetic permeability of the 
atmospheric plasma. At certain photon energies, propagation angles, and matter densities, 
plasma and vacuum contributions cancel each other out, leading to a vacuum "resonance" in 



which both photon mod es become circu 
optics breaks down (see 



Meszaros 



1992 



arly polarized and the approxim ation of geometric 



Lai &: Ho 



2002 



Ho fc Lai 



2003. and the references 



therein). A photon of given energy E, propagating through the density gradient of the NS 
atmosphere, experiences two effects: (1) if it is X mode polarized, it experiences a significant 
enhancement of its absorption opacity in a narrow range of densities about the resonance; 
and (2) partial conversion occurs between X and O states in a manner analogous to the 
MSW effect for neutrinos (ILai &: Hoi 120021 ). The two primary effects of vacuum polarization 
on magnetar spectra are suppression of line features and softening of the hard tail present in 
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previous models (see, e.g.. 



Ho fc Lai 



200ll ) . For details on the ph ysics of vacuum polarization 



and its effect on radiative trans 



(120031), lLai fc Hoi (|2003bJ), and 



er in magnetar atmospheres , see lLai fc Hoi (120021 ) 



Ho fc Lai 



van Adelsberg fc Lail (120061 ). 



It can be shown that photons of energy E, traveling through the vacuum resonance region, 
experience resonant mode conversion with probability Pc = 1 — exp [— tt(E/ E a d) 3 The 
adiabatic energy, E ad , is defined: 



E 



ad 



2.52 



/ B tan0B 1 - (E Bi /Ey 



2/3 



H- 1/3 keV, 



(2) 



where 63 is the angle between the photon propagation direction and the magnetic field, H p 
is the atmosphere scale height, and E^% ~ 0.63(Z/A)B U keV is the ion cyclotron energy, 
with atomic and mass numbers Z and A, respectfully. The parameter fs is a slowly varying 
function of B, whose magnitude is of order unity. At photon energies E <C E^, mode 
conversion is ineffective, while for E > lAE^-, mode conversion is essentially complete. The 
vacuum resonance density occurs at: 



Pv » 0.96 (Z/A) Ef B( A r B l g cm-- 



(3) 



The vacuum resonance phenomenon leav es a unique signature on the emission of polar- 



ized radiation from the NS ( Lai fc Ho 



2003ah . For magnetic field strengths B<7x 10 13 G, 



photons of both modes encounter the vacuum resonance after decoupling from the NS atmo- 
sphere, while for 7x 10 13 G < B < 5 x 10 16 G, the vacuum resonanc e occurs at an atmospheric 
depth between the X and O mode photospheres (lLai fc Holl2003al ). In the former case, pho- 
tons encounter the vacuum resonance while propagating into the magnetosphere. For low 
energy photons with E -C -E^d, mode conversion is ineffective, and the vacuum resonance 
has no effect on the emerging radiation. For high energy photons, mode conversion is effec- 
tive, and results in conversion between the X and O polarization states. Thus, the vacuum 
polarization effect causes the plane of linear polarization to rotate 90° between low and high 
energies. This can occur in models with magnetic fields as high as B = 7 x 10 13 G, though 
there is a dependence on the NS geometry (see §3]). 
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8 van Adelsberg & Perna 
3 EMISSION MODEL 

3.1 Evolution of Photon Stokes Parameters in a Dipole Magnetosphere 

After emerging from the atmosphere, photons propagate through the NS magnetosphere 
before reaching the observer. As discussed above, we ignore scattering and only consider 
magnetic field effects on the polarization mode evolution. 

The solution to Maxwell's equations in magnetized vacuum yields the unit polarizatio n 



vectors for X and O mode photons in the magnetosphere (e.g., Ivan Adelsberg fc Laill2006l ): 

e G = (cos (p B , sin Lp B ), (4) 
e x = (-sm(p B ,cos(p B ), (5) 

where (ps is the azimuthal angle of the magnetic field projected in a plane perpendicular to 
the observer line of sight (see below). Thus, the photon polarization state depends only on 
the magnetic field direction. The difference in the eigenvalues corresponding to the X and 
O polarization states is EAn/(Hc) ~5x lO 3 ^ B\ A (r/R)~ 6 cm -1 , where E x = E/(l keV) 
is the photon energy at the NS surface and r is the distance between the photon and the 
NS. When r>R, the derivative of the angle for a dipole field is approximately dips/ds ~ 
r _1 = 10~ 6 rj^ 1 cm -1 , where r 10 = rj (10 km) and s is the affine parameter defined along the 
photon geodesic. Thus, near the star surface, E An/ '(he) ^> dips/ds, and the polarization 
state evolves adiabatically with the changing direction of the magnetic field. As the photon 
continues to propagate in the magnetosphere, it reaches the polarization limiting radius, r p i, 
where EAn/(hc) = dtps/ds. At distances greater than r p i, the polarization state is fixed, 
and the directions of the mode vectors are constant. Thus, the measured values of the photon 
Stokes parameters are "frozen in" at r p j. Far from the NS, dips/ds ~ l/rj, where ri = c/Q 
is the light-cylinder radius. Therefore, r p i takes the value: 

1/6 



1.5 x 10 3 (ErBlJ^ 1 ) R w km. (6) 



In contrast, if the polarization states of emitted photons were determined near the NS 
surface, the observed signal would be greatly reduced. In this case, addition of the Stokes 
parameters for photons emitted from regions with distinct magnetic field directions would 
tend to cancel. Adiabatic evolution of the photon modes to distances far from the star 
surface, where the magnetic geometry is uniform, leads to significant polarization signals, 



even when emiss ion occurs over an extended region on the star surface ( IHeyl &: Shaviv 



Heyl et al. 



2002 ; 



20031 ). 
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We present calculations for NSs with a pure dipole magnetic field structure. As discussed 
above, realistic magnetar models include twisted magnetic field configurations which contain 
contributions from higher-order multipoles. Nevertheless, we have argued that because the 
polarization signal is fixed far from the star surface, the dipole co mponent of th e field is the 



most important for determining the observed Stokes parameters. iWang &: Lail (120091 ) per- 
formed a detailed analysis of polarization mode evolution in a strongly magnetized vacuum, 
and showed that adiabaticity can be broken near the star surface only when the photon 
traverses a quasi-tangential region (where the photon direction is along the magnetic field). 
They presented a semi-analytic formalism for incorporating this effect into future calcula- 
tions; however, for the polar cap sizes considered in this paper, (3 = 5 — 30°, the correction 



wil 



be small. 



Lai fc Hoi (]2003at ) used the adiabaticity of photon mode evolution in the magneto- 



sphere to develop a simple formalism for calculating the observed photon Stokes param- 



eters. The results of t 



reir ca lculation are accurate fo r NSs with rotatio n periods P > 0.01 s 



( Ivan Adelsberg fc Lai 



20061 ) . Following the work of lLai fc Hoi (l2003al ) , we employ a fixed 
coordinate system, xyz, in which the observer line of sight is along the z axis. The rotation 
axis, fl is defined to be in the xz plane, such that flxz = sinewy, where a# is the angle 
between the rotation and z axes. The NS magnetic dipole moment vector, /a, is inclined at 
an angle ckm relative to the rotation axis. 

We define the rotational phase ip as the azimuthal angle subtended by the magnetic 
dipole vector around the axis of rotation. The phase is taken to be ip — when \x is in the 
xz plane. The angle between [i and z is given by: 

cos0 = cos dji cos Q-m + sina^sinaMCOs^- (7) 

The NS magnetic field, in the near-zone of the star such that r ^ Q/c, can be calculated 
according to the standard dipole formula, B = [3(/x • ?)? — fi] /r 3 . Along the observer line 
of sight, r — z, and the projection of B into the xy plane yields: 

cos (fB = (sin an cos q>m — cos olr sin q>m cosip) /| sin 0|, (8) 
sin(^£ = — sin a M sin0/| sin 0|. (9) 

Far from the NS, the photon geodesic is approximately r « rz. The phase can be written 
■0 = i ) era + r p\/ r u where ip em is the value of the rotational phase when the photon is emitted. 
For magnetars, which rotate slowly, r pl /r; 1, and ip ~ ip em . 
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10 van Adelsberg & Perna 

For rapidly rotating NSs with fi > 100, an accurate calculation of the observed polariza- 
tion requires integration of the transfer equations for the photon Sto kes parameters in the 



magnetosphere, which can result in significant circular polarization (Ivan Adelsberg Lai 



20061 ). Circular polarization can also be generated in cases where th e polarization limitin g 



radius occurs close to either the NS surface or light cylinder radius (jHeyl &: Shavivll2002l ). 
However, for magnetars with adiabatic evolution of the modes, the normalized Stokes pa- 
rameters take the simple form: 

Q E /I E = n| m cos(2^ B ), (10) 
U E /I E = n e i°sin(2^), (11) 

where Ie = + ^f(^em) is the total specific intensity. The Qe, Ue Stokes parameters 

define the plane of polarization in the xyz coordinate sytem. Thus, measurements of the 
photon polarization directly map out the rotation of the magnetic dipole field around the 
NS. 



3.2 Emission from a Finite Spot 



Phase-dependent emissi on from an extended ar ea on a NS surface is compu ted according to 



the m ethod described by 



Pechenick et al 



Perna fc Gottheli 



(11983|), using the generalizations of|] 
( 120081 ). We assume that the radiation emerges from a single spot, centered around the mag- 
netic dipole vector, with opening angle (3. Thus, the angle between the spot center and 
observer line of sight is given by equation ©. We describe points on the NS surface using 
the polar angle 9 and azimuthal angle ip, in spherical polar coordinates. If G = 0, 9 is 



restricted to 9 < /3, otherwise the condition becomes: 
0-/3 <9< + /3 

27T - <p* < if < if*, 

where 

cos (3 — cos 9 cos 9 

cos^ = . . . 

sin B sin 9 

When < < j3, the constraints on 9 and <p become: 

9 < 0*(e,^/3), 

cos/5 = sin 6 sin 9* cos <p + cos 6 cos 9*, 



(12) 
(13) 

(14) 

(15) 
(16) 
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where equation (Tl6]) must be solved numerically. The values of 9 and ip that describe points 
in the spot are restricted to those that yield solutions for co-angles <£>*(/?, 6, 9) or O, ip), 
depending on the condtions. 

Due to general relativistic effects, a photon emitted from colatitude 9 will reach the 
observer if emitted at an angle 5 with respect to the surface normal, where the re l ation 



betw een the two angles is given by the ray-tracing function (IPechenick et al. 



1983 



Page 



1995|): 



9(6) 



R 3 /2R 



du x 



R) \2RJ 



2 2 
U X 



-1/2 



(17) 



In equation (fTTI) . x = sin 5, R s = 2GM/c 2 is the Schwarzchild radius, and M is the NS mass. 
In our calculations, we assume that M = 1.4M and R = 12 km. A simpler, more convenient 
relationship b etwee n the emission angle and colatitude is given by the approximation in 



Beloborodov! (120021 ): 



cos 5 



R s 
~R 



+ 1 



R s 
~R 



cos 9. 



(18) 



Equation (|18|) differs from (fTTI) by less than 1% when i? > 3i? s . The portion of the NS 
surface visible to the observer is given by the set of 9 values that yield solutions for 5 in 
equation (TTBl . 

The total, phase-dependent specific intensities from the spot are then obtained by in- 
tegrating the local emission in each m ode oyer th e observable surface, accounting for the 
gravitational redishift of the radiation ( PagJ 1995 ): 



F j (£oc,V0 



2tt 



dip Ij [9, ip, E^e 



(19) 



where the local specific intensity in each mode Ij (j = X, O) is set to zero outside of the 
boundaries on 9 and ip set above. The energy observed at infinity is given by E^ = Ee A , 
with 



1- 



R s 
~R 



(20) 



The phase-dependent Stokes parameters are then readily calculated using the methods of 



4 RESULTS 

We present results from our magnetar models with fully ionized hydrogen atmospheres, 
magnetic field strengths 7 x 10 13 G - 5 x 10 14 G, and effective temperature T e g = 5 x 10 6 K 
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(a value which is in the typical magnetar range). We explore four representative geometries: 
otR = 90°, aju = 90° (hereafter denoted Gl); an = 90°, «m = 45° (hereafter denoted G2); 
an = 45°, olm — 45° (hereafter denoted G3); and or = 45°, olm = 0° (hereafter denoted 
G4). 

4.1 Atmosphere Model Linear Polarization Fractions 

As discussed above, quiescent magnetar emission spectra are well fit by a blackbody plus 
power-law and do not contain absorption features. Thus, it is difficult to distinguish mag- 
netic atmosphere models from phenomenological fits. However, the linear polarization is 
strongly dependent on the magnetic field and viewing geometry, and, importantly, cannot 
be predicted without a detailed physical model for the surface emission. In this section, we 
describe the dependence of Il em on magnetic field and geometry. 

Figures [I] and [2] show the linear polarization fraction calculated using NS atmosphere 
models at magnetic field strengths B = 5 x 10 14 G (solid curves), B = 10 14 G (dotted 
curves), and B = 7 x 10 13 G (dashed curves). Each panel in the figure corresponds to one 
of the NS geometries described above. The polar cap size is set to (3 = 5°. 

Figure CD shows results for photons with energy E = 0.5 keV. The upper-left panel depicts 
the case of the "orthogonal rotator," in which the magnetic dipole vector rotates in the yz 
plane, intersecting the line of sight at if) — 0. The upper-right and lower-left panels show less 
extreme geometries. In G2, most of the photons are emitted with 5>7r/4, and in G3, the 
dipole vector intersects the line of sight when ■0 = 0, but sweeps out a cone whose base is 
perpendicular to the NS rotation axis. The lower-right panel shows the extreme case when 
the rotational and magnetic axes are aligned. 

There are two key features of Figure [T] that highlight the interplay between the NS 
geometry, magnetic field strength, and vacuum polarization effects to produce the emitted 
polarization. These features include: (1) |II em | is smaller for values of if) close to than for 
values of if> close to 7r/2; (2) at if> ~ 0, the sign of Il em is positive for B = 7 x 10 13 G, and 
negative for stronger magnetic fields. 

Result (1) can be understood by considering the variation in X and O mode opacities 
with angle between the photon propagation and magnetic field directions. For the orthogonal 
rotator, small phases roughly correspond to emission along the observer line of sight and 
magnetic field direction. For emission angles < 5 2tt, the difference between the free- 
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Figure 1. Emitted polarization fraction, III? 1 = C^O — Ix)/(Io + Ix)> as a function of phase, for an extended spot on a NS 
surface. The spot is located at the magnetic pole and has opening angle (3 = 5°. Results are shown for four NS geometries and 
three magnetic field strengths: B = 5 X 10 14 G (solid curve), B = 10 14 G (dotted curve), and B = 7 X 10 13 G (dashed curve). 
The angles an and dim are the polar inclinations of the rotation axis relative to the line of sight, and the magnetic dipole 
vector relative to the rotation axis, respectively. The photon energy is set to E = 0.5 keV. 



free ab sorption opaciti es of the X and O modes is smaller than for angles 5 ~ 7r/4 (see 



52.6 of 



Ho fc Lai 



20011 ). The densities at which X and O mode photons decouple from the 



atmosphere are therefore closer in value, leading to the smaller magnitude of IT 3 ™. Conversely, 
at phases for which ip > tt/2, emission angles for photons reaching the observer are 6 < 7r/4. 
At these angles, the difference between the X and O mode opacities (and hence, decoupling 
depths) is maximal, yielding larger values of |n em |. This trend is also visible in geometries 
G2 and G3 of Figured! where in the G3 case, emission angles 5 ~ tt/4 are realized for phases 
ip w 7T. In the G4 geometry, the emission angle and hence the polarization are independent 
of the phase. 

Result (2) is a manifestation of the vacuum polarization effect. As discussed above, 
for typical angles and photon energies, the emitted specific intensity of X mode photons 
is greater than that of O mode photons. Under these conditions, Il em < 0. However, at 
magnetic field strength B = 7 x 10 13 G, photons with energy E = 0.5 keV encounter the 
vacuum resonance at a density Pv "^Po, where po denotes the O mode decoupling density. 
Thus, if resonant mode conversion occurs, n em will switch sign, corresponding to a rotation of 
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the plane of linear polarization. At emission angles 5 ~ it /A and photon energies E = 0.5 keV, 
the mode conversion probability is much less than unity. However, when the emission angle 
5^0, the mode evolution becomes adiabatic, and Pc — > 1. Therefore, at phases for which 
the emission angle is small, resonant mode conversion occurs and the sign of Il em changes. 
For magnetic fields B > 10 14 G, the vacuum resonance density is intermediate between the 
X and O mode decoupling densities, and no change in the plane of linear polarization occurs. 
Rotation of the plane of linear polarization at B < 7 x 10 13 G occurs in geometries Gl and 
G3 (upper- and lower- left panels), in which %p — > corresponds to 5 —>■ 0. In geometries G2 
and G4, (upper- and lower-right panels), the emission occurs at angles 5 ~ n/4, resulting in 
little mode conversion at E = 0.5 keV. 

As the magnetic field strength is increased, the suppression factor in the X mode free-free 
absorption opacity increases as B\ A . Thus, we expect the emitted polarization fraction to 
increase (roughly) with the strength of the atmosphere magnetic field. When vacuum effects 
are taken into account, the decoupling depth of X mode photons is effectively r educed by 



the in crease in free-free absorption opacity and mode conversion at the resonance (IHo fc Lai 



20031 ). If the X mode photon decoupling density is approximately equal to the vacuum 



resonance density, the trend is preserved (pv oc Bf 4 ). This results in large values of |IT m | at 
B = 5 x 10 14 G relative to those at B = 7 x 10 13 G. 

Figure [2] shows the same results as figure [U, except for photons with E = 2 keV. While 
the general variation in |II em | with phase remains the same as in the E = 0.5 keV case, 
the rotation of the plane of polarization at if) m no longer occurs. This is because, at 
B = 7 x 10 13 G, photons with E = 2 keV encounter the vacuum resonance at a larger 
density, which occurs between the X and O mode photospheres for 5<1, but at a smaller 
density than photons in models with B = 10 14 G. Thus, there is no rotation of the plane of 
linear polarization at the lower field strength, and the magnitude of the polarization fraction 
|n em | is larger for B = 10 14 G. 



4.2 Observed Stokes Parameters 

The observed polarization signals are calculated by multiplying the atmosphere polarization 
fraction by the angular factors in equations ( flOl) and ( TlTl) . Thus, the characteristics 
of intrinsic atmosphere emission are modulated by the projection of the NS magnetic field 
direction into the xy plane. 
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Figure 2. Emitted polarization fraction, III? 1 = C^O — Ix)/(Io + Ix)> as a function of phase, for an extended spot on a NS 
surface. The spot is located at the magnetic pole and has opening angle (3 = 5°. Results are shown for four NS geometries and 
three magnetic field strengths: B = 5 X 10 14 G (solid curve), B = 10 14 G (dotted curve), and B = 7 X 10 13 G (dashed curve). 
The angles an and are the polar inclinations of the rotation axis relative to the line of sight, and the magnetic dipole 
vector relative to the rotation axis, respectively. The photon energy is set to E = 2 keV. 



Figures [3H6] show the observed, normalized photon Stokes parameters Qe/Ie (solid 
curves) and Ue/Ie (dotted curves) at energies E = 0.5 keV and E = 2 keV as a func- 
tion of phase, for NS atmospheres with magnetic field strengths B = 5 x 10 14 G (heavy 
lines) and B = 7 x 10 13 G (light lines). Each panel in the figures corresponds to one of the 
NS geometries described above. In Figures [3] and HI the opening angle of the polar cap is set 
to (3 — 5°, while in figures [5] and E] it is set to (3 = 30°. 

In the case of the orthogonal rotator, sin(2^s) = 0, and cos(2<^b) = —1, yielding 
Qe/Ie = — III", Ue/Ie = 0. Thus, the projection of the magnetic field into the xy plane 
is along the y axis. In general, the projection of the dipole magnetic field into the xy plane 
far from the NS is approximately equal to that of the magnetic dipole vector shifted by tx 
radians (ip — > tps + 7r). It is also important to note that the results are not qualitatively 
sensitive to the polar cap size; the most significant difference between Figures [31 HlandEl M 
are the larger ranges of phase over which the spot is visible. 
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Figure 3. Observed, normalized Stokes parameters Q/I (solid curves), U/I (dotted curves), as a function of phase, for an 
extended spot on a NS surface. The spot has opening angle (3 = 5°. Results are shown for the same geometries as in figure [l] 
for magnetic field strengths B = 5 X 10 14 G (heavy curves) and B = 7 X 10 13 G (light curves). The photon energy is set to 
E = 0.5 keV. 



4.3 Phase- Averaged Polarimetry 

A unique signature of vacuum polarization effects on the evolution of photon modes was first 
reported by lLai fc Hoi (l2003al ). For NS with magnetic fields B < 7 x 10 13 G, the vacuum 



resonance density is less than the decoupling densities of X and O mode photons. Thus, 
photons which have decoupled from the NS atmosphere undergo partial mode conversion at 
the vacuum resonance. At low energies, E <C E^, and Il em is unaffected. However, at high 
energies, E > E a< x, and the sign of n em reverses. Thus, there is a rotation of the plane of 
linear polarization between low and high energy photons. 

In practice, this effect depends on the NS geometry (see §4. ID . For emission angles 5 ~ 0, 
i^ad ~ and no rotation between low and high photon energies is observed. In geometries for 
which the typical emission angle is 5 ~ 7r/4, the effect can be pronounced. Figure [7] shows 
the phase-averaged Stokes parameter Q ave = (27r) _1 / 27r dipQfy) as a function of photon 
energy for several magnetic fields with B = 4 — 50 x 10 13 G. The geometry is case G2 with 
opening angle (3 = 5°. Note that ?7 ave vanishes by symmetry. For B = 4 x 10 13 G (solid 
curve), Q ave changes sign at E ~ 1.5 keV, the energy at which resonant mode conversion 
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Figure 4. Observed, normalized Stokes parameters Q/I (solid curves), U/I (dotted curves), as a function of phase, for an 
extended spot on a NS surface. The spot has opening angle (3 = 5°. Results are shown for the same geometries as in figure [l] 
for magnetic field strengths B = 5 X 10 14 G (heavy curves) and B = 7 X 10 13 G (light curves). The photon energy is set to 
E = 2 keV. 

becomes effective. For B = 7 x 10 13 G, the vacuum resonance density is nearly equal to 
the O mode decoupling density, depending on the NS rotational phase. For most values 
of ijj, p = pv is within the X and O mode photospheres, and no rotation of the plane of 
p olarization occurs. Thus, Q ?V P retains the same sign at low and high energies (c.f., figure 20 



of 



van Adelsberg fc Laill2006l . for a similar result in a different geometry). 



4.4 Polarization, Atmosphere Composition, and NS Equation of State 

Magnetar atmospheres are typically assumed to consist of mostly hydrogen. However, if 
hydroge n is depleted by t hermonuclear burning in the photosphere, they may consist of 



20041 ). For fully ionized helium composition at B = 5 x 10 G, and 



helium (jChang et al. 

T e g = 5 x 10 6 K, we find that the observed Stokes parameters differ from those in the 
hydrogen case by < 10% for 0.1 < E < 10 keV. We therefore do not show results for this 
case. The reason for the similarity is that the O mode decoupling depth in ionized helium 
is ~ 60% of that in hydrogen, while the vacuum resonance density (identified with the X 
mode decoupling depth at large magnetic field strengths) increases by a factor of two. The 
temperature profile evolves much more slowly than the density; in an interval in which the 
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Figure 5. Observed, normalized Stokes parameters Q/I (solid curves), U/I (dotted curves), as a function of phase, for an 
extended spot on a NS surface. The spot has opening angle j3 = 30°. Results are shown for the same geometries as in figure [l] 
for magnetic field strengths B = 5 X 10 14 G (heavy curves) and B = 7 X 10 13 G (light curves). The photon energy is set to 
E = 0.5 keV. 

density increases by several orders of magnitude, the temperature only changes by a factor of 
a few. Thus, at most energies and propagation angles, the relative decoupling temperatures 
of X and O photons are similar to those in the hydrogen case, and we expect to see the same 
qualitative features, with a slightly larger magnitude of the helium polarization fraction. 
Nevertheless, at lower fields, the reduction of the O mode decoupling depth and increase in 
the resonance density supresses the rotation of the plane of polarization between low and 
high energy photons at B « 4 x 10 13 G. For helium, lower magnetic fields, which are unlikely 
to occur in the magnetar regime, are necessary to see this effect. 

We also explored the dependence of the emitted polarization fraction on the equation of 
state of the star. In particular, for the case of B = 5 x 10 14 G and (3 = 5°, we compared 
the polarization fraction for three values of the NS radius, R = 9, 12, 15 km with the NS 
mass fixed at M — 1AM & . For the four geometries defined above, we found that the range 
of ijj for which the spot is visible increases with decreasing radius. This result is expected, 
due to light deflection effects around relativistic stars. However, the difference between the 
emitted polarization fraction for the two extreme radii is |n|p(i? = 15 km) — Il^ n (i? = 
9 km)|/|Il^ n (i? = 15 km)| 1, for most values of ip. Since the phased flux modulation of 
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Figure 6. Observed, normalized Stokes parameters Q/I (solid curves), U/I (dotted curves), as a function of phase, for an 
extended spot on a NS surface. The spot has opening angle j3 = 30°. Results are shown for the same geometries as in figure [l] 
for magnetic field strengths B = 5 X 10 14 G (heavy curves) and B = 7 X 10 13 G (light curves). The photon energy is set to 
E = 2 keV. 

the X and O mode photons is equally affected by changes in the NS radius, the variation 
in the linear polarization fraction with R is associated with two effects that vary with NS 
compactness: (1) the difference between the energy at the star surface and the observed 
energy, and (2) the range of emission angles 6 for which photons reach the observer at a 
given rotational phase ip. The ratio of emitted energies corresponding to the same observed 
energy for stars with R = 9 km and R = 15 km is E 9 /E 15 pa 1.15. Similarly, the difference 
in the range of 5 for photons which reach the observer, at a given 8, is smaller than the 
spot size for small 8, and, for most geometries, is generally smaller than the range over 
which substantial polarization variations occur. Therefore, the polarization signal has no 
significant, direct dependence on the EOS of the NS. 

5 DISCUSSION 

Current studies of NSs rely on spectroscopy to constrain the stars' physical parameters. Fits 
to the phase-averaged X-ray spectrum yield the effective temperature of the star and the 
size of the emitting region, but are insensitive to the viewing geometry. Timing analysis 
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Figure 7. Phase-averaged Stokes parameter Q as a function of photon energy E (keV). 

complements the information obtained from spectral fits by constraining the permitted val- 
ues of an and cum- However, since the amplitude of the modulation of the flux with phase 
is largely dependent on the mass-to-radius ratio of the sta r, there is a degeneracy between 



Perna fc Gottheli 



2008|) 



the inferred viewing geometry and compactness ratio (e.g., 

As shown above, our calculations of the magnetar polarization signal are not sensitive 
to M/R, and can therefore be used to break this degeneracy. Moreover, the sensitivity of 
the Stokes parameters to the magnetic field strength can be used to constrain the magnetar 
magnetic field. In particular, the phase-averaged behavior of the Stokes parameters at low 
and high photon energies in the soft X-ray band can put lower or upper limits on B (in the 
absence or presence of rotation of the plane of linear polarization, respectively). 

Measurements of the magnetic field of a magnetar from the shape of the continuum 
X-ray spectrum alone are extremely difficult to perform. Magnetic atmosphere models are 
characterized by hard tails with respect to blackbody emission (this discrep ancy is larger 



for magnetic fields B < Bq\ see 



Ho &: Lai 



2003 



van Adelsberg &: Lail 120061 ) . Spectral fits 



are generally unable to uniquely discriminate among different B-field strengths based on 
the spectral sha pe alone, especi ally since hard tails are also predicted from scattering in 
the atmosphere (jRea et al.ll2008l ). While no features have been detected in quiescent mag- 
netar thermal spectra, there is one confirmed detection of an abs orption-like feat ure in 
an outburst of the Anomalous X-ray Pulsar 1XRSJ170849-400910 (IRea et all 120031 ). The 
phase-dependent feature occurs at 8.1 keV in the magnetar spectrum, consistent with res- 
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onant cyclotron scattering by electrons at 9 x 10 11 G or ions at 1.6 x 10 15 G. The value of 
the magnetic field inferred from dipole spin-down is B ~ 5 x 10 14 G, however, much higher 
surface fields are possible if the magnetic field is dominated there by higher-order multipoles. 
Nevertheless, an unambiguous identification of the line feature, and hence, magnetic field 
strength, is not possible. Future measurements of Stokes parameters in the soft X-ray band 
may help break this degeneracy. If the rotation of the plane of linear polarization is not 
observed, this is consistent with a the magnetar in the strong magnetic field regime, and the 
identification of the line feature with ion cyclotron resonance becomes more certain. 
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